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Abstract 



In this expository article we describe the asymptotics of certain Fredholm determi- 
nants which provide solutions to the cylindrical Toda equations, and we explain how 
these asymptotics are derived. The connection with Fredholm determinants arising in 
^ the theory of random matrices, and their asymptotics, are also discussed. 

o 

OO ' 1. Introduction 

In the Gaussian unitary ensemble, scaled in the bulk so that the density of eigenvalues 
becomes one, the probability £"(0; s) that an interval of length s contains no eigenvalues is 
equal to the Fredholm determinant of the operator on L 2 (0, s) with kernel 

1 sin7r(x — y) 
it x — y 

This is a result of M. Gaudin ||. This quantity is directly related to the spacings between 
consecutive eigenvalues. The conditional probability density that, given that there is an 
eigenvalue at x, the next larger eigenvalue lies in an infinitesimal neighborhood of x + s is 
equal to —d 2 E(0; s)/ds 2 . This quantity is independent of x since the kernel depends only on 
the difference of the arguments. 

To see how large the gaps between successive eigenvalues might be the asymptotics of 
E(0; s) as s — > oo become of interest. (The asymptotics as s — > are trivial.) A complete 
asymptotic expansion was found by F. J. Dyson ||, and reads 

E(0; s/n) ~ 2 1 / 3 e 3 ^ 1 ) s^ 4 e^ 8 (l + ^ + J + • • •) , 

with computable constants a 1; 02, • • •. The main part of the asymptotics (all but the paren- 
thetical correction terms) were obtained by applying a formal scaling argument to a result of 
the second author |17| on the aymptotics of certain Toeplitz determinants. The asymptotics 
were completed using inverse scattering techniques. 
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The derivation described above was not rigorous. The first rigorous result of this type is 
also due to the second author |1| and was that E(0; s) = e ' s2 / s+ °^. More exactly, 



^log£(0;s) = -~ + O (l). 

This was subsequently extended to a complete asymptotic expansion by Deift, Its and Zhou 
0]. Thus all of the Dyson expansion has now been proved, except for the constant factor. 

Incidentally this constant is expressible in terms of the Barnes G-function which will 
also arise below. This is an entire function satisfying G(l) = 1 and the functional equation 

G(z + l) = T(z)G(z). 

It turns out that 

The kernel (M arising in the Gaussian unitary ensemble has the form 



x-y 



(2) 



with if = sinvrx/v^ and ip(x) = cos'^^x^/7T. Many other kernels of this same form arise in 
random matrix theory and for some there are (heuristically derived) asymptotic expansions 



analogous to the Dyson expansion [|12| , |13|1 . One of the interesting things about the function 
E(0; s) is that it is expressible in terms of a Painleve transcendent, more precisely Py @. 
It was shown by the authors |12], [H| [14| that analogous quantities arising in several other 
matrix ensembles are also expressible in terms of Painleve functions. All the kernels had the 
same general form (Q), and this fact was crucial for the argument. 

The first connection between Fredholm determinants and a Painleve function was es- 
tablished in a paper of B. M. McCoy, T. T. Wu and the first author ||10|| . The function is 
the solution to the cylindrical sinh-Gordon equation (reducible to a special case of the Pm 
equation) 

?"(*) + -?'(*) =8sinh2g (3) 

which satisfies 

q(t) ~ . e~ 4t as t — > oo. 

In the cited work q was shown to be expressible as a certain sum of multiple integrals. This 
sum is in fact equal to a combination of logarithms of Fredholm determinants of two integral 
operators. Soon thereafter an integral representation was found, in terms of this same 
function q, for det (/ — A 2 K 2 ), where K is the integral operator on L 2 (0, oo) with kernel 



X + y 



(4) 



This is also of the form (||). (More precisely it can be brought to that form with the 
change of variable x — > x 1 / 2 .) More recently |3], [TI| it was discovered that there is a direct 
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connection between this function and Fredholm determinants associated with K. This is 
q = log det(I+XK)— log det(J— XK). (With hindsight it can be seen that this representation 
is equivalent to the one found in [fTD| .) It is a generalization of the kernel (Q) which we shall 
be considering here. 

The connection problem for equations such as (§) is to determine the behavior of a singular 
point (in this case t = 0) if the behavior at another singular point (in this case t = oo) is 



known. This connection problem was solved in fllQ| ; the asymptotics will be a special case 
of those we shall describe below. The Dyson expansion, of course, solves the connection 
problem for another Painleve function. 



We shall describe here the results and methods of |L6| which give the asymptotics of the 
Fredholm determinants of a class of kernels containing the above K as a special case. These 
determinants provide solutions to the cylindrical Toda equations 

q'k(t) + \ q' k (t) = 4 (e"- 1 - k e Z. (5) 

If we have a solution q of @ then q k = (— l) k q is a solution to @ which is 2-periodic in the 
sense that g^+2 = <?& for all k. Conversely any 2-periodic solution of @ gives the solutions 
q = (q k - 9fc _i)/2 of (|). 

We introduce a class of n-periodic solutions to the Toda equations generalizing the solu- 
tion stated above for the sinh-Gordon equation. Define a family of kernels indexed by k G Z 
by the formulas 



QJ=£1 



This is clearly ra-periodic in k, and it was shown in [pi that 



q k = log det(J - K k ) - log det(J - K k ^) (6) 

provides a solution to @. In the special case n = 2 there is a single coefficient and we recover 
the solution to (|3|) described above. (The exponential factors in the kernels are different, but 
the determinants are easily shown to be the same.) 

It is the asymptotics as t —>■ of these determinants det(J — K k ) which we shall discuss, 
and of course this will give the asymptotics of the solution q k . The asymptotics take the 
form 

det(I-K k ) ~ bt a , (t->0), 

where the constants a and b are determined by k and the coefficients c u . (We give the 
formulas for them and state our precise assumptions below.) 

Asymptotic formulas of this type in the n = 2 case for det(J — K 2 ) were obtained earlier. 
The constant a was determined in and b in 0. In fact these were obtained for a more 
general class of kernels, where K 2 is replaced by KK 1 and K now denotes the kernel (|j) 
multiplied by a factor {x/y) e . As was already mentioned, the asymptotics of q for the 
n = 2 case were found in ||10|| ; in a method was described to find connection formulas for 



solutions of a class of equations including some of the n = 3 cases of 
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2. Statement of the results 



Since is obtained from Kq by multiplying the coefficients c w by u k we may as well 
assume that k — 0. We do this and write K instead of K . It is also convenient to change 
the constant b so that the asymptotic formula reads 

det(I-K) ~b(-)\ (t->0). 

The constants a and b are expressed in terms of certain zeros cui, • • • , a n of the function 

sin7rs — 7r V] c u {—ujY^ 1 . 

(The exponential is determined by taking | arg(— uS) \ < 7r.) The formulas are 

1 " (n + 1) (2n + 1) 
a = - 2^ a * a ' 

n b1 <ng(^ + i)"- |j| 
n-,=iG(^ + i)- 

Here G is the Barnes G-function mentioned earlier. 

To explain which zeros the «j are, and the assumptions under which this has been proved, 
we introduce a parameter A and the function 



sm7rs 



Denote the zeros of this function by «j(A) (i G Z), where «j(0) = i. Then we take «j = 
for z = 1, • • • ,n. Observe that the «j(A) are analytic functions of A away from those values 
for which the function has multiple zeros, and we choose a path in C from A = to A = 1 
which avoids all such values. Our basic assumption is that there is a path in C from A = 
to A = 1 such that everywhere on the path 

3?ai(A) < 3?a i+ i(A), |3?«j(A) - i\ < 1. 

(We believe that the second condition is redundant.) 

From these formulas we can compute the asymptotics of our solutions (^]) of (|5[). They 
are given by 

q k {t) = Alog(-)+ log B + o{l), 
where for k = 1, • • • , n the constants A and B are given by 

A=2k-fc), 5= n Ls * n 



l<j<k \ n ) k<j<n 1 I n ^ L J 



and for other values of k are given by periodicity. 

In the special case n = 2 this gives the asymptotics for det(J ± XK), with K given by 
(|), and for the corresponding solutions of ([|) for all A except real A satisfying |A| > 1. 
Consequently in this case at least the range of validity is the right one since the asymptotics 



are quite different for these A |T0 
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3. Outline of the derivation 



If R\ denotes the resolvent kernel of K, the kernel of K (I — XK) 1 , then 

d 



Hence 



, ■ logdetfi - XK) = -tr R x . 
dX 

log det( J -K) = - tr R x dX = - / dX R x (x, x) dx, 
Jo Jo Jo 



where the path of integration can be taken to be that in our basic assumption. Now our 
kernel is given by 

£ -t [(1-cj) x+il-iv- 1 ) x- 1 } 

K ( x , y) = Yl c " ; ' 

^ —ujx + y 

and we are looking for the asymptotics as t — > 0. Observe that as t — > 0, 

g-tCi-w- 1 )^- 1 _^ i uniformly for x > 1, e _ ' (1_a;):E — > 1 uniformly for x < 1. 
Therefore if we define kernels K ± (x,y), still acting on L 2 (0, oo), by 

e -t(l-w)a; e -t(l-w- 1 )^- 1 

fsT + (a;,t/) = ■ — , K~(x,y) =^2^ ■ , 

^ —ujx + y ^ —ujx + y 

with corresponding resolvent kernels Rf(x,y) then we might hope that R\(x,x) is close to 
R x (x,x) when x > 1 and close to R\(x, x) when x < 1. This is true, and using operator 
techniques it can be shown that 

J R\{x ) x)dx = J R\{x,x)dx + J R\~{x, x) dx + o(l) 

as t — > 0, uniformly in A. This is heart of the argument. Once we have this approximation 
then everything is plain sailing, since the operators K ± can be transformed into Wiener-Hopf 
operators, and so there are exact integral representations for their resolvent kernels. 

To describe this transformation on K + , say, we observe that we may take t — 1 since the 
variable change x — > x/t replaces K + by the same kernel with t — 1. Assuming this to be 
the case we make the further assumption that each < 0. The reason we can make this 
assumption is that in the present discussion the uj may be arbitrary complex numbers except 
for nonnegative real numbers (when the kernels cease to represent trace class operators) and 
that in the end everything will depend analytically on these uj. So that the formula for the 
resolvent derived with the restriction holds generally. 

Let us define 

A(x, u) = J2^ e- (1 ~ w) x e uxu , B(u, x) = e~ xu . 

These are both kernels of bounded operators from L 2 (0, oo) to L 2 (0, oo) (here is where we 
use 3?cj < 0) and 

K + (x,y)= / A(x,u) B(u,y) du, 
Jo 
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and so K + = AB. Hence (our assumption guarantees that / — XK + is invertible) 

i# = AB(I - XAB)- 1 = A(I - XBA)- 1 B. 
Therefore R\ is known once (J — \BA)~ X is. But the operator BA on L 2 (0, oo) has kernel 

BA(u, v) = -r — ^— — , 

and under the change of variable u — > — 1, t> e 77 — 1 the operator is transformed into 
the operator, also on L 2 (0, oo), with kernel 

£ 



>(£-r?)/2 _ ue -(H)/3' 



This is a Wiener-Hopf operator (its kernel is a function of the difference of the arguments) 
and so its resolvent is expressible in terms of a factorization of its symbol, 1 minus A times 
the Fourier transform of the function 

£ 



,s/2 _ u er m 



oj e~ 

After a change of variable the Fourier transform becomes a Mellin transform and the symbol 
becomes 

suiTrs 

This is how the function (]7]) arises and why everything is expressed in terms of it. 

The symbol is factored into products and quotients of T-functions whose arguments 
involve the zeros «j(A). To see how the Barnes function appears consider an integral of the 
form 



f(s)ds, (8) 



where / has period n. (Integrals like this which involve the a« arise in the computation of 
the constant b.) To evaluate this we integrate 

over the infinite rectangle with vertices | ± ioo and n + | ± zoo. This is equal to @ by the 
periodicity of / and the fact that 

G'(^ + P) G'^ + P) _ r'(* + /3) 
On the other hand this is also equal to 

?G(fT^ res(/ ' a) ' 
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where a runs over the poles of / inside the rectangle. In the integrals that arise these poles 
are exactly the points «i, • • • , a n , and so (§) is equal to 

^gfpf^s(/, Q ,(A)). 
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